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SUMMARY 


This report provides a general analysis of the flight dynamics 
of several surface-to-air and two air-to-air missile configurations. 
The analysis involves three phases: vertical climb, straight and 

level flight, and constant altitude turn. Wind tunnel aerodynamic 
data and full scale missile characteristics are used where available 
unknown data are estimated. For the constant altitude turn phase, 
a three degree of freedom flight simulation is used. Important pa- 
rameters considered in this analysis are the vehicle weight, Mach 
number, heading angle, thrust level, sideslip angle, g loading, and 
time to make the turn. The actual flight path during the turn is 
also determined. Results are presented in graphical form. 



NOMENCLATURE (CONT’D) 


Symbol 
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g 

W 
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X 
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Y 
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01 
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6 

6 
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E 
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Definition 

Weight of propellant grain 
Initial weight of missile 
Weight flow rate of fuel 
Coordinate axis 

Summation of forces on missile along 
X wind axis 

Summation of forces on missile along 
Y wind axis 

Coordinate axis 

Angle of attack 

Sideslip angle 

Elevator deflection 

Rudder deflection 

Atmospheric density 


Heading angle 



NOMENCLATURE 


S3nnbol 
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d 

g 

Isp 

I 

zz 


M 

P 

q 

s 

t 

T 

t. 

V 

w 


bo 


Definition 

Speed of sound 

Drag coefficient 

Zero lift drag coefficient 

Lift coefficient 

Pitching moment coefficient 

Yawing moment coefficient 

Side force coefficient 

Reference length (diameter) 

Acceleration of gravity 

Specific impulse of grain 

Moment inertia of missile about z axis 

Mass of missile 

Mach number 

Atmospheric pressure 

Dynamic pressure 

Body axis yaw rate 

Reference area 

Time 

Thrust 

Time at burnout of grain 
Vehicle speed 
Weight of missile 


V 



I . INTRODUCTION 


In Volume I of this report, aerod 3 mamic coefficients for approxi- 
mately 30 recently declassified missile configurations were compiled 
and tabulated in a consistent format. The purpose of this report, 

Volume II, is to provide a general aerodynamic analysis of several 
missile configurations with maximum utilization of the data in Volume 
I. In this respect, several similar missiles in the surface-to-air and 
air-to-air classes are analyzed, and results are presented in a format 
for performance comparisons . 

Classically, analyses in flight dynamics are divided into three 
major areas. The first is performance, where such items as maximum 
speed, range, ceiling and endurance are considered. These items have 
special forms where tactical missiles are concerned. The second area 
is stability and control. Here maneuverability (control effectiveness) 
and dynamic stability are especially important when coupled with auto- 
matic control systems. Finally, there is aeroelasticity . The high-g 
maneuvers of missiles impose stringent requirements on the aerostructure; 
however, the slender, low aspect ratio configurations in use today are 
relatively Insensitive to classical aeroelastic problems. 

The aerodynamic behavior and performance characteristics of a missile 
are highly dependent on the class of mission for which it was designed. 

In the surface-to-air category five configurations of one missile are 
considered; a wingless version, two similar winged versions, and a 
winged design with aft-tail controls and canard controls. In the air- 
to-air class, two specific designs are examined. 
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There are two approaches for performance anlayses. One Is strictly 
an aerodynamic comparison of different designs with all variables other 
than the aerodynamic coefficients remaining fixed. The other is a 
dynamic comparison incorporating actual variables and initial conditions 
for each configuration, which is the method used herein. Since most of 
the configurations investigated are operational. Information such as 
overall dimensions, weight, and grain size is available. Other variables - 
e.g., specific Impulse and moments-of-inertia - could be reliably estimated. 

As stated previously, the objective of this analysis is to provide 
data on several missiles in certain aerodynamic environments . The 
analysis consists of three phases; two under performance and one under 
stability and control. The two phases in performance Involve straight 
line horizontal flight and straight line vertical cllmib . Constant alti- 
tude horizontal turns due to constant rudder are examined in the stability 
and control phase. 

In each phase, the applicable rigid body equations of motion, in 
the form given by Fogarty and Howe in Ref. 1, are solved numerically on 
the digital computer. The coordinate system used for translational 
displacements is referenced to the flight path and, as pointed out in 
Ref. 1, "...makes much lower accuracy and speed demands on the computer 
than does the body-axis system." For the rotational equations, the 
body axis system is used since again maximum utilization of computer 
time is realized. Assumptions, approximations and details of solutions 
are discussed in subsequent sections. 



II. ANALYSIS 


Perfotmaiice »- Phase I; Horizontal Flight 

The performance phase for horizontal flight consists of two parts. 
In the first part, the thrust required to produce a given steady state 
Mach number Is determined. Initially the vehicle force and moment 
equations were written In the vertical and horizontal planes and were 
solved simultaneously In an Iterative manner. These equations were 


T cos a = qS 

% (M,a) 

(1) 

T sin a + qS 

Cj^(M,a,6g) = W 

(2) 

qSdC^(M,a,6^) 

in Cl 

= 0 

(3) 


Various vehicle weights were considered and the thrust required 
versus Mach number was determined at several altitudes. The results of 
these initial computations indicate that the vehicle weight had little 
effect on the thrust required, as shown In Fig. 1. Consequently, 
Equations (2) and (3) were dropped and the angle of attack was set at 
zero. Equation (1) then reduced to 

T - qS Cjj (M) . (4) 

o 

The second part of the horizontal flight calculations was the 
solution of the simplified (a=0) horizontal equations of motion. The 
purpose of these computations was to provide initial condition inputs 
to Phase m, the constant altitude turn, and to further Illustrate the 
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dependence of the vehicle flight Mach number on the thrust and drag. 
These equations were written as 


X - (T ^ qS (M)/m, 

o 

(5) 

X = ; X dt. 

(6) 

II 

(7) 

W = W - / W dt, 
0 

(8) 


For these equations it was assumed that a=0, W = const = T/Isp and 
that the altitude was constant (sea level) . These equations simulate a 
horizontal launch in which the missile is allowed to accelerate until 
the grain is consumed, at which point the flight is terminated. The 
vehicle flight Mach number is thus determined as a function of the 
vehicle weight or weight of grain consumed. 

In making these computations, the initial vehicle weight and initial 
grain weight had to be estimated. A total vehicle density of 100 Ibs/ft^ 
was assumed. As a further check on the weight estimations, the parameter 
W/S = 750 which was used by Spearman and Fournier in Ref. 2, served as 
a guide. Data given in Ref. 3 was also used as a guideline. As a 
general rule, this number falls within the range of the initial weight 
estimate and the weight of the vehicle after the grain burns out. 

Since portions of these flights were subsonic, subsonic aerody- 
namic data were required. In the majority of the wind tunnel test reports, 
subsonic data were not available, and subsonic drag coefficients had to be 
estimated using the methods of Ref. 4. 
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Performance - Phase II; Vertical Climb 

The vertical climb phase of the analysis consists of vertical launch 
at some constant thrust level. As in the other phases, the Initial 
weight of the vehicle as well as the initial propellant or grain weight 
was required. For each flight the specific impulse of the grain was 
assumed constant at 250 seconds . As a further simplification the weight 
flow rate of propellant was assumed to be 

W = T/Isp (9) 

Obviously, the drag and thrust are quite important in this analysis. 
In order to accurately compute the drag, the drag coefficient was written 
as 


= Cp(M) (10) 

Altitude effects were taken into account through equations defining the 
standard atmospheric density and speed of sound as functions of the 
height Z. Because of the dependence of total aerodynamic forces on both 
Mach number and altitude, maximum thrust levels do not necessarily mean 
maximum altitudes. That is, if the thrust level is high, the missile 
accelerations are high which lead to high drag at low altitudes. On the 
other hand, lower thrust levels may produce lower drags at higher alti- 
tudes. For some of the missile configurations presented in this analysis, 
the thrust levels which produced maximum altitudes were determined. 

Because of the vertical launch from the ground, subsonic drag coef- 
ficient data were also required in this phase. It was also necessary 
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to extrapolate the data to higher Mach numbers than those run In the 
wind tunnel tests . 

The procedure for each flight was to set the thrust at a given level, 
then solve the equations of motion for the subsequent trajectory and 
desired data. The governing equations for vertical flight are 


(W/g) Z = T - qS (M) - W, 

0 

(11) 

q = .5 pa^ M^ 

(12) 

W = W - / W dt , 

(13) 

Z = / Z dt 

(14) 

M = Z/a 

(15) 

Z = h = / Z dt 

(16) 


Flights were made at each of four different thrust levels. The 
thrust levels chosen were designed to cover the range of thrust that 
any one particular missile was assumed capable of producing. During 
each flight the vehicle weight, drag, Mach number and time of flight was 
recorded at specific altitude intervals and appropriate plots were made. 

Stability and Control - Phase III; Constant Altitude Turn 

The third phase of the analysis consisted of a constant altitude 
turn. It is understood that all the missiles considered in this report 
will, in all likelihood, not be required to make a full 180® turn. 
However, the ability of the missile to make a turn is indicative of its 
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abillty to maneuver. Consequently the information derived from a turn 
analysis is very descriptive of the missile aerodynamic efficiency. 

Three different altitudes were selected at which the turns were 
to be made; namely, sea level, 25,000, and 50,000 feet. Thrust levels 
during the turn were constant until the fuel was consumed or until the 
180° turn is completed. Vehicle weight, grain weight, and Mach number 
at the initiation of the turn maneuver was required. In order to be as 
realistic as possible and still stay within the constraints of the over- 
all analysis effort, the weight of the vehicle, corresponding grain 
weight, and Mach number were obtained from the vertical climb phase of 
the analysis. For example, during the vertical climb phase, the missile 
total weight, weight of fuel remaining, and Mach number at 25,000 feet 
were stored for future reference. These three parameters then served 
as initial conditions for the turn at that altitude. This same pro- 
cedure was used for the 50,000-foot altitude turn. 

For the sea level turn, a slightly different approach was taken. 
From Phase I, the missile is allowed to boost horizontally until a Mach 
number of 1.5 is obtained, then the turn procedure is initiated. The 
input weights, vehicle and grain, are obtained from the Phase I analysis 
as previously explained. 

The turn maneuver was begun with the missile centerline on the X- 
coordinate axis. The rudder or control surface to initiate the turn is 
deflected a given amount <5 and the missile then attempts to negotiate 
the turn. The vehicle aerod 3 nramics which dictate the flight path and 
subsequent traj ectory were obtained from the wind txjnnel data. Because 
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of the varying weight, thrust cutoff and c.g. shift, the radius during 
each turn was not a constant. It was assumed that the grain was located 
in the aft portion of the vehicle and consequently the c.g. shifted 
forward as the grain was consumed. This forward c.g. travel and its 
effect on moment of inertia was taken into account in the analysis, since 
stability of the missile is dependent on the location of the c.g. Other 
parameters which had to be estimated were the moments of inertia and 
the damping derivatives C ^ and C The damping derivatives were 
estimated using Lifting Surface and Slender Body Theory techniques. 

In the initial solutions of the governing equations, the angle of 
attack, a, was included as a variable and was computed along with other 
parameters. However, it was observed that under the constraints of 
the constant turn maneuver, a had a negligible effect on the resultant 
flight path or other pertinent variables. Consequently, a was eliminated 
in the remainder of the computations. Initial efforts to use linear 
aerodynamics led to Incorrect results; consequently, values of Cy vs. 6 
and Cji vs. 3 were input in tabular form to Include nonlinear effects. 

Double interpolation algorithms were set up in order to determine 
accurately the side force coefficients and yawing moment coefficients. 
Drag coefficients, as in Phases I and II, were extrapolated to M = 1.0 
and M = 6.0 and were also linearly interpolated from input tables. 

As previously mentioned, the thrust level for any one flight was 
held constant at one of four preselected values. At the beginning of 
each turn, the burnout time for the grain was determined from 


t 


bo 


= W /W . 

g 


(17) 
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It was then assumed that tailoff (final burning of grain) occurred 
over a two-second period and the thrust linearly decayed to zero during 
this time. 

Since the turn was at a constant altitude, motion was restricted 
to the x-y plane. The missile was assumed to "slide" through the turn 
(no roll) with the fins or cruciform wings in the "plus" configuration. 

With these restrictions and assumptions, the resultant equations 
of motion were written: 


X = T cos a cos 8 
w 

+ qS(-Cjj cos B + sin B) , 

(18) 

Y = - T cos a sin 
w 

B + qS(Cp sin B + cos B), 

(19) 

^ - 


(20) 

dt . 


(21) 

ii = R^/cos a , 


(22) 

= / iji dt , 


(23) 

<■ 

II 

< 


(24) 

V = / V dt , 


(25) 

B = Y /(mV) - 
w 


(26) 

B = / B dt , 


(27) 


Cy = CY(e.M) 


(28) 
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S - S W <29) 

O 

= C^(3,M) . (30) 

These equations (18) through (30) were solved digitally and the desired 
parameters plotted. The output from the simultaneous solution of 
these equations Includes (1) X versus Y, (2) g-level during turn, 

(3) during turn, (4) $ during turn, (5) Mach number during turn, 

(6) weight of vehicle during turn, and (7) time required to make the 
turn. 



Ill . RESULTS 


The table below Is a general summary of the missiles which were 
analyzed as described In the preceding sections. 


TABLE 



PHASE I 

PHASE II 

6=10“ 

PHASE III 
6-20“ 

6=40“ 

TM X-2774 

X 

X 

X 



TM X-1751 

X 

X 

X 

X 

X 

TM X-1025 

X 

X 

X 

X 

X 

TM X-3070 

X 


X 

X 


TM X-846 

X 


X 

X 


TM X-2780/A 

X 

X 

X 

X 

X 

TM X-2780/C 

X 

X 

X 

X 

X 


For each missile configuration plots are presented which depict the 
results from each of the three "missions.” On each graph is an identi- 
fication number which is indicative of certain parameters as described 
below. Example: 


1005 TM X-2774 



Rudder ^ Thrust in Missile identification 

deflection 1000 lbs. or Report Number 

in degrees 


For this example the rudder deflection is 10 degrees, the thrust is 
5000 lbs., and the missile is taken from NASA Report Number TM X-2274. 
Data on each graph is presented for three different altitudes; sea 
level, 25,000 ft., and 50,000 ft. 

- 11 - 
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In addition to the identification number the figures are numbered 
relative to the respective "mission, " The Roman numeral is respective 
of the mission as outlined below. 

I - Horizontal Flight 
II - Vertical Flight 
III - Constant Altitude Turn 

Other numbers in the figure number designation are in consecutive order. 

As a typical example of the data presented for each missile, the 
missile designated TM X-2774 is discussed in detail. From the hori- 
zontal flight phase Fig. 2-1 depicts the terminal Mach number versus 
thrust levels for various altitudes. For example, if the thrust were 
10,000 pounds, the maximum Mach number the missile could obtain at 
sea level would be about M=3.3 but at 20,000 feet altitude the missile 
could reach a Mach number of M=5.5. 

Continuing with Fig. ''-II, the time required to climb vertically 
to various altitudes is presented. Note that on this particular missile 
the maximum altitudes are very high - in excess of 100 miles - and con- 
sequently the trajectories or data presented : ' these high altitude 
regimes may be somewhat in error. In Fig. 4-II, the Mach number - 
Altitude trace is presented. Note that maximum Mach numbers range from 
5.5, at 5000 lbs, thrust, to 6.5 at 20,000 lbs. 

The last figure in the Phase II flight set is the vehicle Weight 
versus Altitude. Fig. 5-II depicts the weight of the vehicle during 
the vertical boost. From this graph the initial weight of the vehicle 
may be obtained as well as the final burnout weight. 
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The third phase flight, constant altitude turn, is presented in 
the next series of curves. In Fig 6-III through 9-III, the flight 
time and vehicle weight are presented as functions of the downrange 
distance, Y. The initial weight of the vehicle at the constant turn 
altitudes are obtained from Fig. 5-II. The first of the figures. 

Fig. 6-III, is for a thrust level of 5000 lbs. while the second, third 
and fourth are for 10,000, 15,000, and 20,000 lbs. thrust, respectively. 

In order to get a clear picture of the missile performance, all 
the figures in this phase should be considered since they are all 
interrelated. The Mach number during the flight and "heading angle," 
ip, are presented in Figs. lO-III through 13-111. Note that the initial 
Mach number for each altitude is different. The initial values for 
the vehicle flight Mach number is taken from the vertical flight phase. 
Fig. 4-III. Note that for this particular missile and thrust = 5000 lbs. 
at altitude of 50,000 feet, the grain bums out before the 180° turn is 
completed. However, for the other two altitudes, fuel still remains 
after the turn is completed. Increasing the thrust level as shown on 
Figs. ll-III through 13-III produces somewhat different results. The 
next set of curves. Figs. 14-IV through 17-III, present the g loading 
and sideslip angle, 6, as functions of the downrange distance, y. Note 
that at the beginning of the turn maneuver, the missile oscillates 
about some nominal sideslip angle but then damps out and maintains a 
near constant 6. 

The final set in the Phase III series is the most descriptive for 
the constant altitude turn. In Figs. 18-III through 21-III the flight 
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path, X versus Y, for each thrust level is presented. Small replicas 
of the missile shape are drawn on the flight path trajectory at discrete 
locations at its correct heading and sideslip angle. The flight for 
each of the three different altitudes are shown on these figures . 

The same general figures as discussed above are presented for the 
other missiles as outlined in the Table. There are a few exceptions; 
for example, there is no Ph^se I or Phase II flights for TM X-3070 and 
TM X-846 since these are air launched air-to-air missiles. The inter- 
pretation of the results of these computations are left to the reader. 
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Fig. 65-III. Heading Angle and Mach No. vs. Downrange Distance, Y 
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Fig. 85-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 94-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 96-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 97-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 100--III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 116-III, Heading Angle and Mach No. vs, Downrange Distance, Y. 
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Fig. 135-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 137-III. Heading Angle and Mach No. vs. Downrange Distance, Y 
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Fig. 145-III. Constant Altitude Flight Path, X vs. Y. 
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Fig. 152-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 163-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 


178 



/ 




Fig. 164-III. Flight Tine and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 165-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 167-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 203-III. Heading Angle and Mach No. vs. Downrange Distance, Y 
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Fig. 206-III. Number of G*s and Sideslip Angle vs. Downrange Distance, Y. 
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Fig. 210-III. Constant Altitude Flight Path, X vs. Y. 
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Fig. 214-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 216-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 217-III, Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 218-Iir. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig* 220'-III. Heading Angle and >Iach No. vs. Downrange Distance, Y. 
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Fig. 221-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 222-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y. 
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Fig. 223-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y 
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Fig. 224-I1I. Number of G's and Sideslip Angle vs. Downrange Distance, Y 
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Fig. 225-III. Number of G’s and Sideslip Angle vs. Downrange Distance, Y 
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Fig. 226-IH. Constant Altitude Flight Path, X vs. Y. 
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Fig. 228-III. Constant Altitude Flight Path, X vs. Y 
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Fig. 236-II1. Heading Angle and Mach No, vs. Dovmrange Distance, Y. 
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Fig. 240-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y 
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Fig. 242-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y. 
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Fig. 243-III. Constant Altitude Flight Path, X vs. Y. 
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Fig. 24A-III. Constant Altitude Flight Path, X vs. Y. 
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Fig. 248-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 249-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y 
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Fig. 251-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 252-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 254-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 255-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y. 
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Fig. 256-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y. 
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Fig. 257-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y, 
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Fig. 258-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y, 
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Fig. 260-III. Constant Altitude Flight Path, X vs. Y. 
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Fig. 264-III, Flight Time and Vehicle Weight vs. Downrange Distance, Y 
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Fig. 265-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y 
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Fig. 268-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 269-III. Heading Angle and Mach No. vs. Downrange Distance, Y 
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Fig. 271-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 272-III. Number of G’s and Sideslip Angle vs. Downrange Distance, Y. 
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Fig. 274-III. Number of G's and Sideslip Angle vs. Downrange Distance, Y. 
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Fig. 275-III. . Number of G*s and Sideslip Angle vs. Downrange Distance, Y 




Fig. 276-III. Constant Altitude Flight Path, X vs. Y. 
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Fig. 281-III, Flight Time and Vehicle Weight vs, Downrange Distance, Y. 
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Fig. 282-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 


297 



WEIGHT-LBS m10‘ TIME-SEC 

J3.00 50.00 100.00 150.00 J).00 20.00 40.00 60.00 


CD 


,2015TMX-846 



Fig. 283-III. Flight Time and Vehicle Weight vs. Downrange Distance, Y. 
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Fig. 284-III. Heading Angle and Mach No. vs. Dovmrange Distance, Y 
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Fig. 287-III. Heading Angle and Mach No. vs. Downrange Distance, Y. 
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Fig. 291-III. Number of G's and Sideslip Angle vs. Dovmrange Distance, Y. 
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Fig. 295-III. Constant Altitude Flight Path, X vs. Y. 
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